Charge fluctuations in a final state with QGP 
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Abstract 

Charge fluctuations as a possible signal of quark - gluon plasma (QGP) were 
recently suggested. A short summary of comments presented on this subject is given 
and supplemented by a discussion of the coexistence of pions produced "directly" 
and through a QGP phase. Such a coexistence may obscure the expected plasma 
signal similarly to the effects considered in the comments mentioned above. 
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It was suggested recently [|], that the "primordial frozen fluctuations" may be used 
to signal the formation of a quark-gluon plasma (QGP) in the early stage of high energy 
interaction of heavy ions. Similar considerations led to more detailed predictions for 
measurable quantities P, p|. In particular, Jeon and Koch || argued that the event-by- 
event fluctuations of charge (or, equivalently, of the ratio of positive to negative pions) in 
a restricted rapidity range may be used as a tool to signal the possible formation of quark 
- gluon plasma (QGP). They compared the ratio of charge dispersion squared to charged 
multiplicity in two models: the "pion gas" and the QGP and conclude that they differ by 
a factor of five. Similar result were found for the dispersion of "positive-to-negative" ratio 
for pions. The dramatic difference may be easily understood as the reflection of small 
quark charges as compared to hadrons and of zero gluon charges. 

Jeon and Koch concluded that the strong decrease of charge fluctuations in the QGP 
as compared to a pion gas should be seen as an "unmistakable signal of QGP formation 
from 'Day-1' measurements" at RHIC. They added some caveat about resonances and 
other correlation effects which may reduce the fluctuations in the "pion phase", but these 
were claimed to be minor corrections. 

After appearance of this paper we indicated || that the values of the ratio of charge 
dispersion squared to charged multiplicity R quoted in Ref.|| are unrealistic. In any 
reasonable model of standard (non-QGP) hadron production these values depend strongly 
on the rapidity bin width and already for the values of a few units they are as low as 
those quoted for QGP. This was demonstrated on the example of the JETSET/PYTHIA 
generator ||. This dependence stems mainly from the effect of global charge conservation, 
which was discussed simultaneously elsewhere 0. 

This prompted the authors of Ref.|| to introduce corrections in the journal version 
of their paper M and to write an extended paper, in which the role of the corrections 
for global charge conservation is discussed in more detail ||. It was suggested that after 
dividing the measured ratio R by a correction factor 1 — p, one recovers the originally 
suggested values for "QGP" and "non-QGP" cases. | Thus the authors still regard their 
analysis as providing a clear signal of possible QGP formation. 

This claim was questioned in a few papers. It was pointed out that the predictions of 
string models of "non-QGP" hadroproduction for the R ratio differ quite strongly from 
those of thermodynamical models [[10]. This makes the possible distinction of "QGP case" 
rather questionable. Also the results from different MC generators were analyzed in detail 



lTJ. It was argued that the possible independence of R on the impact parameter may 
suggest that it measures a QCD scale rather than QGP effects. Another effect neglected 
in Refs. M, |Sfl is the final state rescattering, which may reduce significantly the difference 



between the "QGP" and "non-QGP" values of R g| . 

In this note we consider one more obvious effect, which may impede the observation of 
a QGP signal: the coexistence of hadrons from two sources in the same event. Indeed, it 
seems rather unrealistic to assume that all the pions observed in the event come from the 
hadronization of QGP. Even for the most central heavy ion collisions some hadrons are 
likely to be produced in the collisions of "peripheral" nucleons, which may be described 
by string- or thermodynamical models. In other words, it seems unlikely that the "QGP 



•fBy p we denote here the probability that a produced pion falls into the considered rapidity bin. 
Obviously, p is an increasing function of bin width. 
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bubble" covers all the interaction volume. 

Let us consider a toy model in which there is a common multiplicity distribution 
P(N) for all pions produced in full phase space. Let us assume that there are binomial 
probability distributions for this pion to be a "direct" positive or negative pion or an 
"effective quark" which later hadronizes in a pion. We assign charges ±1/3 to these 
quarks (which is a rough average for realistic quarks and gluons) and impose global charge 
conservation for produced particles (neglecting thus fragmentation of initial nucleons, 
which does not affect the spectra in the central rapidity region). We find then 

P(N + , N_, K + , K_) = ]T P(N)B q (N, N + )B q (N, N_)5 N;N++ n.+k ++ kJn + -n.,(k.-k + )/3 

N 

^ ( (1) 

where N + , N_ denote "direct" positive and negative pions, K + , K_ positive and nega- 
tive quarks (which hadronize later into pions) and B q (N,M) is the binomial probability 
distribution given by 

W tf »° Ml/M)! ' M "-'»" < 2 > 

One should add that we assume implicitly that each "effective quark" hadronizes into 
one charged hadron. In the old-fashioned recombination model, in which there is one 
charged hadron per 2/3 of qq pairs the suppression effect from QGP phase would be 
much weaker. 

Now let us consider the rapidity interval of the width A. Assume that for each type 
of particles the probability to fall into this interval is given by the same function of bin 
width A (denoted by p, as before) and that the probability to find a given number of each 
type of particles is also given by a binomial formula. Then 

p(n + , n_, k + , fc_) = £ P(N + , AL, K + , K_)B P (N + , n + )B p (N_,n-)B p (K + , h + )B p (K_, fc_). 

(3) 

This allows us to calculate the generating function of charge distribution in the rapidity 
interval 

G%(z) = £^P(Q) = Y,z n+ ~ n - +k+/3 - k - /:i P(n + ,n-,k + ,k^. (4) 
Q 

We get 

G%(z) = G[(l -p + z l ^p) l ' 2 (l -p + z'^pf/ 2 * (5) 
q(l -p + zp)(l ~P + z-^ 3 p) q(l-p + p/z)(l-p + z 1 / 3 p) 

* [ q+ (l-p + Z 1 /^) 2 ){ q (l-p + Z-V3 p )2 )l 

where Q[y\ is the generating function of the multiplicity distribution in the entire phase 
space 

g[y} = J2y N P(N). (6) 

N 

The derivatives of the generating function taken at 1 are the distribution moments. 
In particular 

^-\ V=1 = Y,NP(N)=<N>. (7) 

dy N 
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For the moments of the charge distribution we find then 



<Q> = ^=1 = 0, (8) 

<Q 2 > = <Q(Q-l)>=^\ z =i=<N>^±p(l-p). (9) 
Thus the ratio of charge dispersion squared to average charged multiplicity is 

fl(A) = -^ = i5i±i[l-p(A)]. (10) 

Obviously, for q — 1/2 (which corresponds to all pions produced directly) we recover 
the standard result Dq / < n c h > = 1 — p, and for q = (when all pions come from 
plasma) we get a suppression factor of 1/9 []. For the realistic values of 2q around 0.5 
(which corresponds to a similar number of pions from both sources) the suppression factor 
is above 1/2 instead of 1/9. 

Our toy model represents quite well the broad class of models in which the final state is 
a superposition of pions from two independent sources. In fact, the results hardly change 
if we impose charge conservation separately for both sources. Therefore we may conclude 
that the presence of pions which do not come from plasma brings the suppression factor 
much closer to 1 than expected for QGP. If the fraction r of pions comes " directly" , and 
the suppression factor for plasma is s, we get approximately 

D 2 

R(A) = 9_ = r r + (i_ r ) s ][i_ p (A)]. (11) 

< "./- > 

As one sees, we get the same 1 — p factor for both contributions from "non-QGP" 
and "QGP", as suggested in Ref. 0. This is, however, because we assumed the same 
dependence of p on A for "QGP" and "non-QGP" contribution, which is by no means 
obvious. In fact, it is more likely that the pions from QGP are more peaked at mid- 
rapidity. Then q and r increase with A, which may cancel the decrease due to the 1 — p 
factor. This suggests that the change of the dependence of R on A may be a similarly 
(in) significant signal of the onset of QGP as the drop in the value of R. 

Summarizing, we have shown that the coexistence of the "direct" pions and pions 
coming from plasma reduces the suppression effect expected from "pure QGP". Adding 
this reduction to other effects considered before (resonances and string effects reducing 
the R value for "non-QGP", final state rescattering enhancing the value for QGP etc.) 
one may conclude that it is unlikely to find a strong suppression signal when the threshold 
for QGP formation is passed. 
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1 This very small value is a direct consequence of our assumptions for "effective quark" charges. In a 
model taking into account " true" quarks and gluons this factor is usually 1/5 to 1/4 || . 
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